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ABSTRACT 

We study three processes that eject hypervelocity (> lO^kms"^) stars from the 
Galactic center: (i) close encounters of two single stars; (ii) tidal breakup of bi- 
nary stars by the central black hole, as originally proposed by Hills; and (iii) three- 
body interactions between a star and a binary black hole (BBH) . Mechanism (i) 
expels hypervelocity stars to the solar radius at a negligible rate, ~ 10~^^yr~^. 
Mechanism (ii) expels hypervelocity stars at a rate ~ 10~^{r]/0.1) yr~^, where rj is 
the fraction of stars in binaries with semimajor axis Ob ^ 0.3 AU. For solar-mass 
stars, the corresponding number of hypervelocity stars within the solar radius 
-Ro = 8kpc is ~ 60(77/0. l)(ab/0.1 AU)^/^. For mechanism (iii), Sgr A* is assumed 
to be one component of a BBH. We constrain the allowed parameter space (semi- 
major axis, mass ratio) of the BBH. In the allowed region (for example, semimajor 
axis of 0.5 X 10"'^ pc and mass ratio of 0.01), the rate of ejecting hypervelocity 
stars can be as large as ~ 10~^ yr~^ and the expected number of hypervelocity 
stars within the solar radius can be as large as ~ 10^. Hypervelocity stars may 
be detectable by the next generation of large-scale optical surveys. 

Subject headings: black hole physics - Galaxy: center - stellar dynamics 



1. Introduction 

The Milky Way, as well as most nearby elliptical galaxies and disk galaxies with spheroids, is 
believed to house a massive black hole (BH) at its center (e.g. Schodel et al. 2002; Genzel et al. 
2003; Schodel et al. 2003; Eisenhauer et al. 2003; Ghez et al. 2003a,b,c; Gebhardt et al. 2003). 
The deep gravitational potential well close to the BH enables relativistic phenomena, such 
as the relativistic motion of jets associated with some active galactic nuclei (e.g. Begelman, 
Blandford & Rees 1984). Hills (1988) first pointed out that stars may also be expelled from 
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the vicinity of the BH with high velocities. The main features of Hills' argument can be 
explained simply. Consider a star with specific energy E = ^v"^ + $(r), where v and $(r) are 
the velocity and specific potential of the star at position r. Close to the BH, <l>(r) ~ —GM,/r 
where G is the gravitational constant and M, is the BH mass. Now suppose that the star 
approaches close enough to the BH so that |$(r)| ^ In this case, the star will have 
a velocity v = ^1\E - $(r)] ~ ^j2GM,/r = 2.9 x 10^kms-\M./10^MQy/^{lmpc/ry/^, 
where mpc = 10~^pc will be often used as a length unit in this paper. If the star then 
suffers a velocity change 5v <^ v (e.g. due to interactions with surrounding stars or the 
tidal disruption of a close binary) and the increase of the specific energy of the star SE — 
^{v + Svy — \v'^ ~ v5v is much larger than \E\, the star will escape from the BH with velocity 
roughly given by 

^2v5v ~ 1.5 X lO^kms-^(M./lO^M0)^/^(lmpc/r)^/^(5u/4OOkms-^)^/^ (1) 

The discovery of such hypervelocity stars — stars with velocity exceeding 10'^ km — would 
provide strong evidence for the existence of a massive central BH in our Galaxy — if any 
is still needed — based solely on stellar kinematics in the outer parts of the Galaxy. The 
properties of hypervelocity stars would also illuminate the nature of the innermost regions 
of galaxies, such as the stellar kinematics, the age and metallicity distribution, the number 
and mass of the BHs, etc. 

In this paper, we will study the expected properties of such hypervelocity stars. We 
consider three possible scattering processes. 



1. Gravitational encounters of single stars. The gravitational interaction of two stars of 
masses rrii and rrij leads to a velocity change of star m^, 

2Gm, 



Sv — 



2Gmi 



[2GK + m,-)yd'/' 

where Yij is the impact parameter and Wij is the initial relative velocity of the two 
stars (Binney & Tremaine 1987). Only rare close encounters (impact parameters of a 
few solar radii) can produce velocity changes large enough for a star to escape from 
the BH at high velocity. 

2. Tidal breakup of a binary star. We denote the component masses of the binary as mi 
and 7712. A binary with semimajor axis Ob will be tidally disrupted by the massive BH 
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when it passes closer to the BH than the Roche or tidal radius ~ ab[M,/ (7711+777.2)]''^'^^ ^ 
10 AU(ab/0.1 AU)[M,/10^(mi+7n2)]^^^- During the tidal breakup, each star will receive 
a velocity change of order its orbital velocity relative to the center of mass of the binary 
star, that is, the velocity change of mi is about (Hills 1988; Gould & Quillen 2003) 



'G{mi + 7772) / rn2 



Ob \777i+7772 

er,^3- r^^V'V^V'V^V (3) 



5v 



mi + 7772 / V 1 / V 

If one star gains energy during this process, the other will lose energy and become more 
tightly bound to the BH. Comparing equations (2) and (3), we see that the semimajor 
axis of the binary star Cb in equation (3) plays the same role as the impact parameter 
Yij in equation (2). 

3. Ejection by a binary black hole (BBH). This is actually a special case of mechanism (1), 
in which the perturbing star 777^ is one component of a BBH. The existence of BBHs 
in some galactic centers is a necessary consequence of hierarchical galaxy formation 
(Begelman, Blandford & Rees 1980; Yu 2002; Milosavljevic & Merritt 2002; Volonteri, 
Haardt & Madau 2002; Komossa 2003)^. We denote the component masses of the BBH 
as Ml and M2 (Mi > M2). If the semimajor axis a, of the BBH is smaller than a 
certain value ah(= GM2/4(T^, where is the one-dimensional velocity dispersion of 
the galactic center, see Quinlan 1996, cq. 29), most of the low angular momentum stars 
that pass through its vicinity will be expelled with an energy gain after one or several 
close encounters with the BBH, and the average velocity increase of these stars is given 
by: 

VMi + M2y \IQ^Mq) \ a, j ' ^ ^ 

where the force per unit mass from the BH M2 is F ~ GM2/al and the interaction time 
St ~ [G{Mi + M2)/al]-^/'^ (Quinlan 1996). Note that dv in equation (4) is independent 
of the stellar mass, while dv in equations (2) and (3) depends on the mass of the ejected 
star; otherwise the three equations have similar form. 

This paper is organized as follows. In §2, we review the data on the inner structure of 
the Galaxy, such as the mass of the central BH, stellar density, etc. In §3, we study the rates 



^We do not consider the possibility of three or more massive BHs in galactic centers, although this 
configuration may also occur in some galaxies. 
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of ejecting hypervelocity stars due to the three processes discussed above. A discussion and 
conclusions are given in §4. 

2. The Galactic center 

Measurements of stellar radial velocities and proper motions at radii r > 10 mpc, the discov- 
ery of variable X-ray emission around Sagittarius A*, and the tracing of the Keplerian orbits 
of stars at radii r ~ 0.1-1 mpc have provided strong evidence for the existence of a massive 
BH in the center of the Galaxy (Genzel et al. 2000; Baganoff et al. 2001; Schodel et al. 2002; 
Gcnzel ct al. 2003; Schodel et al. 2003; Eisenhauer et al. 2003; Ghez et al. 2003a,b,c). We 
expect that the distribution of old stars should be approximately spherical and isotropic, 
since the two-body relaxation time is short, and this expectation is consistent with the ob- 
servations (Gcnzel et al. 2000). We may describe the mass distribution near the center of 
the Galaxy by the combination of a point mass of 3.5 x 10^ Mq (obtained by averaging the 
values (4.0 ± 0.3) x 10^ Mq in Ghez et al. 2003c and (2.9 ± 0.2) x 10^ Mq in Schodel et al. 
2003, Fig. 11), plus the visible stellar system with mass density 

p{r) = Ar~" exp{—r^/rl). (5) 

Here a — 1.8, the characteristic radius of the bulge is — 1.9 kpc, and the normalization 
A can be obtained by setting p{r = 0.16 pc) = 3.6 x IO'^MqPC"^ (Dehnen & Binney 1998; 
Schodel et al. 2003). The density distribution in equation (5) has a cusp at small radii, 
p oc r~" as r — s> 0, consistent with theoretical models of a stellar system with a central BH 
and all stars having the same mass {a = 1.75, e.g. Bahcall & Wolf 1976). At the smallest 
radii, r < 10" = 0.4 pc, infrared observations suggest that a decreases to 1.4 (Genzel et al. 
2003), somewhat below the theoretical models. 

Throughout this paper we shall assume that the distance to the Galactic center is 
Rq = 8 kpc (Eisenhauer et al. 2003). If not otherwise stated, we assume for simplicity 
that the stars in the Galactic center have the solar mass and radius. The effect of the 
generalization to a distribution of masses and radii will be discussed briefly in §4. 

We define the distribution function (DP) /(x, v) of a stellar system so that /(x, v) d'^xd'^v 
is the number of stars within a phase-space volume d^xd^v of (x, v). By Jeans's theorem, 
the DP in a spherical and isotropic stellar system depends on (x, v) only through the spe- 
cific energy E. The DP f{E) in the Galactic center, which will be used to calculate the 
ejection rates of hypervelocity stars below, can be obtained from equation (5) and the Ed- 
dington formula (eq. 4-140a in Binney & Tremaine 1987). The potential in the stellar system 
$(r) = $,(r) + $*(r) is contributed by both the central BH ($,) and the bulge stars ($*); 
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however, at small radii r < 1 pc the gravitational potential contributed by stars can be 
ignored. In this region we have $ = = — GM./r, p{r) oc r~", and the stellar DF is 

f{E) = C\Er'/\ (6) 

where 

C = 2-3/V-3/=m±^ . a(r) . J^K, (7) 

T{a — ^) V r 

and is the stellar mass. Note that C is independent of radius since the radial variations 
of p(r) and cr(r) cancel. 



3. Rates of ejecting stars with high velocities 

In this section, we study the rates of ejecting stars with high velocities from the Galactic 
center. We will consider in turn the three ejection mechanisms described in §1. 



3.1. Close encounters between two single steirs 

Two-body gravitational encounters between single stars can eject one of the two stars (Henon 
1969; Lin & Tremaine 1980). Consider a test star with mass rrii and velocity Vj that en- 
counters a star with velocity Wj. The test star will experience a velocity change from Vj to 
Vj + (see eq. 2 for |ejj|). If the test star is moving through a background of stars with DF 
/(r, v), the probability that this star will suffer an encounter with a background star during 
time dt, rj(r, Vi)dt, is given by: 

r,(r, v,)di ^dt j dV,- /(r, v,)«;,,- j dFj j d*,,-, (8) 

where Yij is the impact parameter, Wij = |vj — \j\ is the relative velocity and \E'jj is the 
angle between the (vj, v^) plane and the (vj — Vj,eij) plane (Henon 1960). If the velocity 
distribution of field stars is isotropic, we have 

r,(r, v,)di ^dtj dvj ^f{r,vj) J dwij J dV^^ J d*,,-, (9) 

where Vj = The probability that the test star is expelled with a velocity higher than 

Vq at infinity can be obtained by restricting this integral to the encounters that change the 
value of its velocity from Vi = |vi| to a value larger than v^, i.e., 

\Vi + eij\ > Ve- (10) 
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Here ^v^ + $(r = oo) = ^v^ + $(r) = Eq, and we shall assume that $(r = oo) = 0. We 
call vo the ejection speed and £'0 the ejection energy. Thus, the rate of ejecting stars with 
ejection speed higher than vq due to two-body encounters in a unit volume around position 
r, 7^(r, Vo), is given by: 

7^(r,^;o) = J d\^ f{r,Yi)Ti{r,Yi) 

= STT^y J dt. y d«;,,y" drj J d*,,-/(r,v,)/(r,v,.)^i^.-<-, (H) 

where the integration variables should satisfy the inequality (10). The total rate of ejecting 
stars with ejection speed higher than vq is given by 



dt 



47r / n{r,vo)r^dr. (12) 







The ejection rate can be obtained analytically in some cases. Ignoring possible stellar 
collisions (when the impact parameter Yij is small enough) and treating stars as point masses 
with zero radius, equation (11) can be simplified as follows (Henon 1969; Lin & Tremaine 
1980): 

7^(r,.o) = ^TT^GV fdEfiE) f ^E' f{E') ^^' f' ^l'^^"'" . (13) 

By ignoring the gravitational potential contributed by stars (e.g. in the region close to the 
BH) and applying equation (6) in equations (12) and (13), the total ejection rate of stars as 
a function of ejection speed can be found to be (Lin & Tremaine 1980): 

diVej(^;.) _ 2VV/2G'^M.V /£.V""'/'r2(l + «)r(|-2a) ^^^^ 



dt 3 (t9 V^V (3-q;)(2-q;) ' 

Although the ejection rate in this equation is integrated over radius, we have expressed it in 
terms of the velocity dispersion a{r) and the density p(r), given by equation (7) — these can 
be evaluated at any radius since their radial dependences cancel out in equation (14). Note 
that in this formula the ejection speed v, represents the velocity after escaping the gravity 
of the central BH but before climbing through the bulge, and the ejection energy E, — ^v^. 
Using a — 1.8, BH mass M. = 3.5 X 10^ Mo, and the stellar mass density distribution 
p{r) = Ar~" (see eq. 5 and the subsequent text) in equation (14), we obtain the ejection 
rate shown in Figure 1 as the dotted line at the top of the figure. 

In reality, stars have non-zero radii, and encounters of two stars with impact parameter 
Yij < Ycoii = {Ri + Rj) [1 -I- 2G{mi + rrij) / {Ri + Rj)WijY^'^ will lead to stellar coUisions, where 
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Ri and Rj are the radii of the two stars rrii and rrij. Such coUisions have two principal effects 
on the calculation of the ejection rate: (i) Encounters with Yij < Ycon should not be included 
in equation (11). (ii) Stars cannot diffuse to energies E < —Gm^/R^ without colliding first, 
where i?* is the stellar radius (Bahcall & Wolf 1976). Thus, we cut off the stellar distribution, 
so that f{E) = for E < —Gm^/R^. After making these modifications, the ejection rate 
is hard to obtain analytically so we resort to numerical calculations, shown by the dashed 
line in Figure 1. As seen from Figure 1, the ejection rate is significantly decreased by stellar 
collisions, especially at speeds > 300 km s~^ that approach or exceed the escape speeds from 
the stellar surface. 

These results do not include the effects of the gravitational potential of the bulge, disk 
and halo. To include the effects of the bulge, we use equation (5) to compute the gravitational 
potential from the bulge; we add the potential from the BH and use Eddington's formula to 
compute the bulge DF f{E)] we then truncate the DF for E — $*(r = 0) < —Gm^/R*, and 
insert the result into the formula (13) to obtain the ejection rate. The result is shown by the 
bold solid line in Figure 1; this is much lower than the dashed line at low speeds, because 
low-velocity stars cannot climb the potential to escape the stellar bulge. The difference 
between the bold solid line and the dashed line is negligible at ejection speeds higher than 
1000 kms^^ 

These results do not include the potential of the Galactic disk and halo. The gravita- 
tional potential difference between Rq = 8kpc and the outer parts of the bulge, rb — 2kpc, 
can be roughly estimated by 



where is the circular speed, assumed independent of radius. The velocity distribution of 
the stars that reach the solar radius Rq is given by 



The thin solid line in Figure 1 shows the ejection rate as a function of Vh, which is obtained 
from the bold solid line using equations (15) and (16). 

The ejection rate of stars at speeds > lOOkms"^ is about (1-2) xlO~^yr~^ at the outer 
edge of the bulge and about 2 x 10~^ yr~^ at the solar radius; the ejection rate of stars at 
speeds > 500 km s^^ is only 10"^ yr^^ at the edge of the bulge and a factor of three lower at 
the solar radius; and the ejection rate at speeds > 1000 km s^^ is only 10^^^ yr^^ at either 
location. We conclude that the expected rate of formation of hypervelocity stars in the 
Galaxy due to two-body close encounters of single stars is negligible. 



A$ = y^Mn(i?o/rb) ~ |(370kms-^)2(yc/220kms-^)^ 



(15) 




(16) 
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Finally, to illustrate the importance of the BH, we show the ejection rate from a stellar 
system with the same density distribution as that in the Galactic center but without a central 
BH (the dot-dashed line in Figure 1). The existence of a central BH dramatically increases 
the rate of ejecting hypervelocity stars. 



3.2. Tidal breakup of a binary star by the central BH 

It is plausible that binary stars are formed as commonly near the Galactic center as they are 
in the solar neighborhood. Close encounters between a close binary and a massive BH may 
break up the binary through an exchange collision, in which one binary component becomes 
bound to the BH and the other is ejected at high velocity (up to several thousand kms~^; 
Hills 1988; Gould & Quillen 2003). 

We shall focus on binaries in the scmimajor axis range 0.01 AU < Ob < 0.3 AU; 0.01 AU 
is nearly 2Rq and thus is a lower limit to the binary separation for solar-type stars, and the 
choice of an upper limit of 0.3 AU will be explained later. The binary period distribution 
observed locally by Duquennoy & Mayor (1991) implies that roughly 8% of binaries belong 
to this period range; given that the fraction of stars in binaries is ~ 0.5 we expect that the 
binary fraction in the semimajor axis range 0.01 AU < % < 0.3 AU is about 0.04. This 
estimate is quite uncertain, both because the estimate derived for the solar neighborhood is 
uncertain and because we do not know whether the binary distribution is the same at the 
Galactic center as it is in the solar neighborhood. 

We must first check that these binaries are not disrupted by encounters with single stars. 
The disruption time is (Binney & Tremaine 1987) 

where p is the density of stars (eq. 5), Ob is the scmimajor axis of the binary, a is the velocity 
dispersion. In A is the Coulomb logarithm, and we have assumed that the binary components 
and the field stars have the same mass. In the region r < 1 pc where the gravitational 
potential is dominated by the BH, we may take a ~ {GM,/ry/'^ and A ~ {M,/m^){ab/r), 
so that 

4 X W^jr f r /O.l AU~ 



ln[1.7(ab/0.lAU)(lpc/r)] VlpcJ V «b J ' ^^^^ 
This formula is derived for soft binaries, for which A > 1; a similar result holds for hard 
binaries, except that In A is replaced by unity, and tjis represents the timescale in which the 
semimajor axis shrinks by a factor of two due to encounters. We conclude that binaries with 
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Fig. 1. — Estimated rates of ejecting stars from the Galactic center by two-body encounters 
of single stars. Here v, is the velocity after escaping from the BH, vq is the velocity after 
escaping from the bulge, and Vh is the velocity after reaching the solar radius. The dotted 
line (a function of v,) is obtained by treating stars as point masses with zero stellar radius 
and ignoring the gravitational potential contributed by bulge stars (see eq. 14). The dashed 
line (a function of v,) is obtained by assuming that all stars have the solar radius. The bold 
solid line (a function of ■;;o) is the ejection rate including the effects of both non-zero stellar 
radii and the gravitational potential of the bulge. The thin solid line (a function of Vh) is 
the rate at which ejected stars reach the solar radius (see eqs. 15 and 16). The current of 
hypervelocity stars with velocity > 700kms~^ past the solar radius is neghgible (less than 
1 per Hubble time). The dot-dashed hne is the ejection rate including the effects of both 
non-zero stellar radii and the gravitational potential of the bulge, but assuming that there 
is no central BH. The existence of the BH increases the ejection rate, especially at high 
velocities. 
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semimajor axis Cb = 0.3 AU can survive for the age of the Galaxy outside r ~ 0.8 pc, while 
binaries with Cb = 0.01 AU can survive outside r ~ 0.05 pc. 

We now review Hills' (1988) three-body simulations of the encounter between a binary 
and a massive BH. The center of mass of the binary is initially on an unbound orbit relative 
to the BH, with velocity Voo- The probabihty of an exchange coUision is a decreasing function 
of the dimensionless closest- approach parameter 

here wc have assumed for simplicity that mi = 1712 = ni^, and -Rmin is the closest approach of 
the binary to the BH (the periapsis of the orbit of the center of mass of the binary around the 
BH). The parameter R!^i^ is essentially the distance of closest approach in units of the Roche 
radius of the binary. Hills finds that exchange coUisions occur in nearly 80% of encounters 
with ~ 0.3 and in 50% of encounters with 7?^;^ ~ 1. Let us define an ejection speed 
parameter by (eqs. 1 and 3) 



where v, is the rms velocity at infinity of the ejected star. Then v'^ is mainly a function of 
and roughly independent of the BH mass and the semimajor axis and masses of the 
binary components. The velocity v', varies between 130 and 160kms~^ for R'.^^^ — 0-1, with 
a peak near i?^in — 0-3, and is about 130kms~^ at i?min = 1 (Hills 1988). 

These results are independent of the initial relative velocity Voo of the binary relative to 
the BH. This assumption is correct so long as this is much less than the relative velocity at 
-Rmin, which in turn requires that^ 



Ob / V"^ 



For the Galaxy this limit corresponds to v^o ^ 1-4 x 10^kms~^(0.1 AU/ab)"*^^^, which is 
satisfied by a large margin. Similarly, Hills' results for the probability of an exchange collision 
should be vahd for orbits that are bound to the BH so long as the orbits are very eccentric, 
that is, i?min <^ o whcrc a is the semimajor axis of the orbit around the BH. 



^The dependence of the exchange cross-section on Uoo is mapped out in more extensive three-body inte- 
grations by Hills (1991), although Hills' empirical summary of the results is based on the assumption that 
the critical velocity scales as mI^^ , not mI^^ as in equation (21). 
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For numerical estimates we assume that both binary components have the solar mass. 
According to equation (19), for a binary to have an exchange probability of 50% or more in 
an encounter with a 3.5 x 10^ Mq BH, it must pass within about -Rmm — 150ab of the BH; 
that is, within -Rmin = 1-5 or 45 AU for = 0.01 and 0.3 AU respectively. The rms ejection 
speed at infinity v, in this case is about 5 x 10^ or 0.9 x lO^kms"^ for ab=0.01 or 0.3 AU 
respectively. These ejection speeds for Ob < 0.3 AU are high enough to escape the Galactic 
bulge and even the Galactic halo, which is why we restrict our attention to this range of 
binary semimajor axis. 

The ejection rate depends on the flux of binary stars passing close to the central BH, 
specifically, within a periapsis distance i?min — 150ab. The region in the {J,E) =(specific 
angular momentum, specific energy) phase space where a binary star has dimensionless 
closest approach parameter R'^^^^ ^ 1 is called the "loss cone" and is given by: 

< Jl{E, Ob) = 2Rl,^ [E - *(i?„,in)] 2GM.R^^ {\E\ < GM,/R^^). (22) 

[Note that the concept of the "loss cone" is used in other contexts, for example, see Magorrian 
& Tremaine (1999) for the tidal disruption of single stars by a massive BH or Yu (2002) and 
§3.3 for the expulsion of low-angular momentum stars by a BBH. The loss cone for disruption 
of single stars is usually much smaller than the loss cone for disruption of binary stars.] At 
first, the ejection rate is determined by the rate of depletion of the initial population of 
binary stars in the loss cone, and the rate per unit energy is given by (see eq. 19 in Yu 2002): 

F'-\E, ab) ~ 47rV(^)^ic(^, «b), t< P{E), (23) 

where P{E) is the radial period of an orbit with energy E and zero angular momentum, and t] 
is the number fraction of binary stars. As the binary stars in the loss cone are disrupted, new 
binaries arc scattered into the loss cone by two-body relaxation. Eventually, the scattering 
into the loss cone reaches a steady state. The steady-state diffusion rate of binary stars 
into the loss cone per unit energy is obtained from the Fokker-Planck equation (eq. 24 in 
Magorrian & Tremaine 1999): 

^ rgFf-»(E,Gb)/ln[GM./(4|E|i?^g], if g « - ln( J^/ J,^) 



F'^\E,a^), ifg»-ln(j2/j2)' 
where 

q = P{E)^{E)Jl{E)/Jl{E, ab). (25) 

Here Ji is the orbit- averaged diffusion coefficient and Jl{E) is the specific angular momentum 
of a circular orbit at energy E. The quantity \qJic is the mean-square change in the scalar 
angular momentum for a nearly radial orbit in one orbital period; \qJl = P{E){{AJf) = 
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P{E){r'^{AvtY), where (•) denotes an orbit-averaged diffusion coefficient and Vt is the tan- 
gential velocity (see details in Appendix B of Magorrian & Tremaine 1999). The region in 
which q > —ln( Jj^/J^) is known as the "pinhole" or "full loss cone" regime, as opposed to 
the "diffusion" or "empty loss cone" regime. Generally the loss cone is full at large radii and 
empty at small radii, we shall define the critical radius rcrit to be the apoapsis of orbits at 
the transition between these two regimes. The critical radius rcrit ~ 600 pc for Ob = 0.3 AU 
and ~ 20 pc for = 0.01 AU. 

We can use the model for the Galactic center described in §2 to estimate the ejection 
rate due to disruption of binary stars. As usual we assume that the DF is isotropic in velocity 
space, and that it vanishes for energy E' — $^,(r = 0) ~ ii^ < —G Mq/ Rq (without this cutoff 
the ejection rate for a full loss cone diverges). If the binary loss cone is full, the ejection rate 
is given by equation (23): 



J F^'^\E) dE 

1-0 

~ 8Ti^r]GM,R^inC / \E\'^-^/ME 

---^--(i^)(a^)(¥)(al^)"^°' (-) 

which is approximately consistent with the estimate ~ 10~^(?7/0.1)(ab/0.1 AU) in Hills 
(1988), based on much less certain observational parameters. However, this calculation, 
like Hills', neglects the fact that the binaries initially in the loss cone are rapidly depleted. 
For a more accurate ejection rate we must use equation (24), which yields 

Uei = J F,,{E) dE ~ 1.5 X 10-^(7^/0.1) yr-i (27) 

if Ob = 0.3 AU. The ejection rate is quite insensitive to the binary semimajor axis a^, 
decreasing to only rZej = 0.9 x 10-^(7^/0. 1) yr~^ if Ob = 0.01 AU. The apoapsides of most of 
the stars contributing to the total diffusion rate into the loss cone are in the range 2-3 pc 
for Ob = 0.01-0.3 AU; at these radii the binary disruption rate due to encounters with field 
stars (eq. 18) is small. Because the loss cone is not full. Hills' simple estimate of the ejection 
rate is too large by three orders of magnitude at Cb = 0.1 AU. 

The time required for the ejected stars to travel a distance D is by: 

D , _6 f D 1.6 X 103kms-i\ 

tj, ~ _ = 5xl0Vr^T 28 

V, \8kpc V, / 



5 X 10^ yr 



D 1.3xl02kms 



-1 



8 kpc v'. 



-13- 



M. 'Mo; 1,0.1 AU'mJ ' ^^^> 

where v', is the ejection speed parameter defined in equation (20). Prom equations (27) and 
(29), the number of hypervelocity stars inside the sphere with radius D is given by: 



1/3 



T] \ f D \ / ab \ 1/2 / M© 



In this process, one component of the binary star is ejected with an energy gain, and the 
other loses energy and becomes more tightly bound to the BH. The apoapsis distance of the 
bound star is approximately GM,/vl ~ 1 x 10^ AU(M,/3.5 x 1O6M0)(1.6 x 10=^ kms-Vu.)^. 
The eccentricity of the bound star is usually high, and the corresponding orbital period is 
~ 8yr(M,/3.5 x 1O^M0)(1.6 x 10^ kms""^/i',)^. Recent high-resolution infrared observations 
of the region around the compact radio source (and presumed BH) SgrA* have revealed 
massive young stars on high-eccentricity orbits with periods of decades or less (Schodel et 
al. 2002, 2003; Eisenhauer et al. 2003; Ghez et al. 2003a,b,c). Gould & Quillen (2003) have 
suggested that these unusual orbits may have been produced by Hills' mechanism; that is, 
these stars may be one component of a tidally disrupted binary star. According to equation 
(27), we expect that the number of young (10^ yr old) stars that are the remnants of tidally 
disrupted binaries should be ~ 100(77/0.1) at any time. Of course, the number of high-mass 
remnants would be smaller, by an amount that depends on the shape of the initial mass 
function. 



3.3. Interactions between a single stetr and a meissive BBH 

BBHs are likely to exist in some galactic centers (Begelman, Blandford & Rees 1980; Yu 
2002), and perhaps in the Galactic center as well. In fact, Hansen & Milosavljevic (2003) 
have recently argued that the young stars observed at 0.1-1 mpc distances from SgrA* may 
have been shorn from a disk or cluster surrounding the smaller component of a BBH, by 
tidal forces from the larger component. A massive BBH steadily loses orbital energy to the 
surrounding stellar population. When the two BHs are orbiting independently in the galaxy, 
and even after they form a bound binary system, the energy loss is mainly through dynamical 
friction. Once the BBH becomes sufficiently tightly bound that its orbital velocity exceeds 
the velocity dispersion of the stars, dynamical friction becomes less and less effective, and 
three-body interactions with low-angular momentum stars passing through the BBH become 
the dominant energy-loss mechanism. When the BBH becomes "hard" , at a semimajor axis 
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(Quinlan 1996) 

4(7,2 Vo.l/ V3-5 X lO^Moy V (^c / 

where 

t/ = M2/(Mi + M2) (32) 

is the mass ratio of the small BH to the BBH, most of the low-angular momentum stars 
that pass near a BBH will eventually be expelled with an energy gain after one or several 
encounters with the BBH. In such expulsions, the average ejection speed of the stars is 

^ / 3.2GMM ^ 2 2 X 10^ km s- r ^ V ( ) (33) 

(see eq. 17 in Yu 2002). The average relative energy change AS/S per ejection is independent 
of the BBH orbital energy S. As the semimajor axis of the BBH continues to decrease, 
eventually gravitational radiation takes over as the dominant energy loss mechanism. 

We now ask what are the observational constraints on the mass ratio u and semimajor 
axis a, of a possible BBH in the Galactic center. We shall assume that the total mass 
Ml + M2 = 3.5 X 1O^M0, that the BBH is on a circular orbit, and that Sgr A* is the larger 
component Mi of the BBH. (We shall argue below that Sgr A* is unlikely to be the small 
component of the BBH.) For reference, the orbital period of the BBH is 

P = 27ra^/^ ^ / 3.5xlO^Me \ 



There are several constraints on the properties of the BBH: 



• A census of stars within ~ 10" ^ 0.4 pc of the Galactic center reveals that the peak 
of the stellar surface density agrees with the position of Sgr A* within ~ 0'.'2 ~ 8mpc 
(Genzel et al. 2003). Taking the center of mass of the BBH as the center of the stellar 
distribution, and ignoring the small possibility that the two BHs are aligned along the 
hne of sight to us, we have the following constraint on the distance between Sgr A* and 
the center of mass of the BBH: 



i/a. < 0'.'2 ~ 8mpc. (35) 

In panel (b) of Figure 2, the region below the bold long-dashed line gives the parameter 
space represented by inequality (35). 



- 15 - 



• The eccentric Keplerian orbits of stars at radii r ~ 0.1-1 mpc (Schodel et al. 2002, 2003; 
Ghez et al. 2003a,b,c) will eventually provide strong constraints on the properties of a 
BBH at the Galactic center. However, a binary with semimajor axis small compared to 
the periapsis distance of the stellar orbit, or large compared to the apoapsis distance, 
will be poorly constrained because the orbit will be nearly a Keplerian ellipse, with one 
focus at either the center of mass or one component of the BBH, respectively. Moreover, 
so far these stars have only been observed for a fraction of an orbital period, so the 
sensitivity to deviations from a Keplerian orbit is still relatively small (but growing 
rapidly). Without carrying out proper simulations, we guess that the agreement of 
the observations to date with orbits around a single point mass constrains the mass 
fraction to be z/ < 0.2 if a, is in the range ~ 0.1-1 mpc. According to this constraint, 
the BBH is unlikely to be in the box enclosed by the dot-long-dashed line near the 
right boundary of panel (b) of Figure 2. 

• Observations of the proper motion of Sgr A* also constrain the parameter space (a,, v) 
of a putative BBH in the Galactic center. VLBA observations from 1995 to 2000 show 
that the dominant term in the proper motion of Sgr A* with respect to extragalactic 
radio sources comes from the motion of the Sun around the Galactic center. The 
rms residual from uniform motion is A ~ 0.5 mas = 0.02 mpc, and the upper limit 
of the peculiar motion of Sgr A* perpendicular to the Galactic plane is v± ~ 8kms~^ 
(Backer & Sramek 1999; Reid et al. 1999, 2003). We do not know the motion of the 
Sun in the Galactic plane with comparable accuracy, but unless the peculiar velocity 
of Sgr A* happens to be closely aligned with the Galactic plane, it is unhkely to be 
larger than \/2v± ~ Vmax — 12kms~^. If the orbital period of the BBH Pbbh (eq. 
34) is short compared to the observation interval T = 5 yr, then the amplitude of the 
motion of Sgr A* relative to the center of mass of the BBH must be smaller than the 
rms observational error, that is, 

i^a, < /A, (36) 

where / is of order unity. In the opposite limit, if Pbbh ^ T then the velocity of 
Sgr A* relative to the center of mass must be less than f max, so 

^ /t' max Pbbh/ (27r). (37) 

We set / ~ 2 and simply choose the transition between inequalities (36) and (37) so 
that they are continuous. In panel (b) of Figure 2, the region on the left side of the 
bold solid line gives the parameter space allowed by inequalities (36) and (37) with 
A = 0.5 mas and Wmax = 12kms~^ (Reid et al. 2003). 

The rate of ejecting stars by interactions with a hard BBH depends on the flux of low 
angular momentum stars passing through the vicinity of the BBH. Just as in the case of the 
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loss cone of binary stars discussed in §3.2, here the region in the (specific energy, specific 
angular momentum) phase space in which a star can pass within ~ a, of the BBH is called 
the loss cone, and may be found by replacing -Rmin with the semimajor axis of the BBH 
in equation (22) (see Yu 2002). Initially, the loss cone of the BBH is full. As stars are 
removed from the loss cone by close encounters with the BBH, new stars refill the loss cone 
by two-body relaxation (or by non-axisymmetric gravitational forces from the surrounding 
stellar system, if it is non-spherical; see Yu 2002). Eventually the system reaches a steady 
state controlled by the balance between the ejection rate of stars and the rate at which 
stars refill the loss cone. If the rms angular momentum transfered to or from the stars per 
orbital period is larger than Jic (see eq. 22), then the stars will refill the loss cone as fast 
as it is depleted and the loss cone remains full; otherwise, the stars will slowly diffuse into 
the loss cone and the loss cone remains nearly empty. The refilling rate due to two-body 
relaxation can be found by solving the steady-state Fokker-Planck equation (e.g. eq. 17 in 
Magorrian & Tremaine 1999 or see Yu 2002). If the BBH is hard and has a total mass 
3.5 X 10® Mq, our calculations show that the rate of removing stars from the loss cone is 
^BBH — 2yr~^(a,/lmpc) when the loss cone is full and ngf^ ~(l-3) x 10"'^ yr~^ when the 
loss cone is empty; the superscript "diff" refers to the fact that the refilling of the loss cone 
is a diffusive process. The rate u^bbb. is insensitive to the BBH mass ratio and the rate 
^BBH insensitive to both a, (only through the logarithmic term in eq. 24) and z/, so long 
as the BBH is hard and the mass of the smaller component of the BBH is significantly larger 
than the stellar mass (Quinlan's 1996 simulations only cover the range u > 0.004 but his 
finding that the hardening is largely independent of u should extend to more extreme mass 
ratios). The hardening timescale of the BBH is then given by (see eq. 33): 



thia,) 



a, 
a. 



_ 6 X lO^r 1^ - (38) 

3m,<^ff ^ V 3.5 X IO6M0 / V / V nf'^ /' ^ ^ 



''BBH V'-'-'-' ^ ^^10/ y / y "BBH 

which is shown as a function of a, and u in Figure 2 (see the solid horizontal lines in panels 
c and d, which represent 4, 6, 8 x 10^ yr from top to bottom). Assuming that the BBH orbit 
is circular, its gravitational radiation timescale is given by (Peters 1964): 

5 c^at r A .f^s f «. Y {^-5x10^ Mq)^ f 0.1 



a. 



~ 5.4xl0Vi' 



64 G^MiM2{Mi + M2) ■ \lmpc J (Mi + M2)2Mi V 

(39) 

which is also shown in panels (c) and (d) of Figure 2 (the long-dashed lines, which rep- 
resent 10^, 10^, 10^ yr from top to bottom). A similar diagram to panel (c) is shown in 
Hansen & Milosavljevic (2003). The BBH loses energy mainly by interacting with stars 
when igr(a., J^) > th{a,) and mainly by gravitational radiation when t^{a,, u) < th{a,). 

If we do not live at a special time, the lifetime of the BBH orbit, min(ih,igr), should 
exceed u x 10^° yr (i.e. at most, the central BH should accrete a fraction u of its mass in 
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a fraction u of its age). This constraint is shown as the bold long-dashed hne in panels (c) 
and (d) of Figure 2. We may combine this with the constraints from the position of Sgr A* 
(inequality 35), the near-Keplerian orbits of stars close to Sgr A*, and the proper motion of 
Sgr A* (inequalities 36 37), to conclude that the BBH should not be located in the dotted 
region in panel (d) of Figure 2. Note that wc have assumed that Mi coincides with Sgr A* 
and is the larger component of the BBH, so that cannot exceed 0.5 (eq. 32). However, 
if we relax this constraint, that is, if v is allowed to exceed 0.5, we may extend the above 
constraints on the BBH parameter space and find that the forbidden region in Figure 2(d) 
excludes the range 0.5 < u < 1. Thus, Sgr A* cannot be the smaller component of the BBH. 

The bold dot-short-dashed line in Figure 2 (panels a and d) represents the semimajor axis 
at which the BBH becomes hard, given by equation (31). The average ejection velocity v, as a 
function of a, and i/ is shown by short-dashed lines in these panels (these represent v, — 400- 
2000 km s~^ with interval Sv, — 200kms~^ from top to bottom, with v, — 1000 km s~^ shown 
in bold). Some short-dashed lines with v» > 1000 km s~^ cross the allowed region of panel 
(d), and hence there are allowed parameters for a BBH that could eject hypervelocity stars 
that reach the solar radius. 

Thus, for example, consider a BBH with u — 0.01 and a, = 0.5 mpc. This has orbital 
period Pbbh — 0.6 yr, gravitational radiation timcscalc tgr ^ 4 x 10^ yr, and distance of 
Sgr A* from the center of mass of 0.1 mas, consistent with the observational constraints 
we have described. The time for a hypervelocity star to travel D = 8kpc is tjy ~ 8 x 
10^yr(Z}/8kpc)(10"^kms^"'^/f,) (cq. 28). Since ti) ^ min(th,^gr), wc can ignore the BBH 
orbital evolution, so the number of stars with ejection speed higher than lO^kms^^ within 
the sphere D < Skpc can be roughly estimated to be u'^bu'^d — 10^. 

4. Discussion and conclusions 

We have studied stellar dynamical processes that eject hypervelocity (> lO'^kms^^) stars 
from the Galactic center. We consider three mechanisms of ejecting stars: close encounters 
of two single stars, tidal breakup of binary stars, and three-body interactions between a star 
and a BBH. 

The rate of ejection from close encounters of single stars is strongly affected by the finite 
size of stars: Figure 1 shows that this rate is 10^ times smaller for solar-type stars than for 
point masses with the same radial distribution. For solar-type stars, the rate of ejection with 
speed > lO^kms"^ at Rq — Skpc is only 10~^^yr~^, which is too low to be detectable. If 
there is a distribution of stellar masses, the ejection rate is hkely to be larger for the fighter 
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Fig. 2. — In each of the four panels, the horizontal axis represents the mass ratio u = 
M2/{Mi + M2) of a hypothetical BBH in the Galactic center and the vertical axis represents 
its semimajor axis a,. SgrA* is taken to be the more massive component of the BBH. 
The eccentricity of the BBH orbit is assumed to be zero. Panel (a): the bold dot-short- 
dashed line gives the semimajor axis of the BBH when it becomes hard (see eq. 31, with 
(Tc = lOOkms"^). The short-dashed lines represent the average ejection speeds of the stars 
from the BBH (eq. 33), v, = 400-2000 km s^^ from top to bottom with interval 200 km s^^, 
with V, = 1000 km s^^ in bold. Panel (b): the region above the bold long-dashed line is the 
BBH parameter space that is excluded by the location of SgrA* close to the center of the 
stellar cusp (see inequality 35). The rectangular region enclosed by the bold dot-long-dashed 
hne is excluded by the Keplerian orbits of stars at radii r ~ 0.1-1 mpc. The region on the 
right side of the bold solid line is excluded by the proper motion of Sgr A*, except for unusual 
orientations in which its motion is almost along the line of sight (inequalities 36 and 37). 
The kink at a, ~ 5 mpc represents an approximate boundary between the cases in which 
the period of the BBH is long or short compared to the present duration of proper-motion 
observations. Panel (c): the long-dashed lines give the gravitational radiation timescale of 
the BBH: 10^, 10^, 10^ yr from top to bottom (eq. 39). The solid lines give the hardening 
time th. 4, 6, 8 X 10^ yr from top to bottom (eq. 38). The region above the bold long-dashed 
hne represents the constraint that min(th,tgr) > x 10^° yr, so that the larger BH accretes 
at most a fraction 1/ of its mass in a fraction 1/ of its age (i.e. the present accretion rate of BH 
companions is not unusually large compared to the average over the lifetime of the Galaxy). 
As in panel (a), the bold dot-short-dashed line gives the semimajor axis of the BBH when 
it becomes hard. Panel (d): combination of the curves in panels (a)-(c). The shaded region 
is excluded. See §3.3. 
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stars (Henon 1969), but these also have lower luminosity and thus are harder to detect. 

The rate of ejection from breakup of binary stars (Hills 1988) depends on the orbital 
parameters of the binaries, such as the semimajor axis and component masses. We estimate 
that the rate of ejecting hypervelocity stars is 10~^(77/0.1) yr"-*^ (eq. 27), where rj is the binary 
fraction and the result is only weakly dependent on semimajor axis for Ob ^ 0.3 AU. This 
rate is also not strongly dependent on the stellar mass distribution. With the assumption 
that the binary components have the same mass m^, the number of hypervelocity stars within 
the solar radius is ~ 100 for rj = 0.1 and Ob = 0.1 AU (eq. 30). 

SgrA* may be one component of a BBH (Hansen & Milosavljevic 2003). We discuss 
the constraints on the parameter space (mass ratio and semimajor axis) of the BBH using a 
variety of observational and theoretical arguments. The ejection speed of stars due to three- 
body interactions with the BBH depends on the BBH semimajor axis and mass ratio, but 
not on stellar mass. In part of the allowed parameter space, near a, = 0.5 mpc and u = 0.01, 
the average stellar ejection speed can exceed lO^kms"^ (see panel (d) of Fig. 2), and the 
rate of ejecting hypervelocity stars is ~ 10~^yr~^. In this case the expected number of stars 
with ejection speeds higher than lO^kms""^ within the solar radius r < 8kpc is ~ 10^. 

Nearby high- velocity halo stars bound to the Galaxy cannot have velocities that exceed 
the local escape speed. Recent estimates of the escape speed are 550 ± lOOkms"^ (Leonard 
& Tremaine 1990), 610 ± 120 km s"^ (Kochanek 1996), and 600 km s"^ (Wilkinson & Evans 
1999), well below our cutoff at lO^kms"^, so it is unhkely that bound halo stars can appear 
to be hypervelocity stars. Rare Local Group interlopers, perhaps from M31, could have 
velocities exceeding the escape speed. However, hypervelocity stars ejected from the Galactic 
center have an additional distinctive kinematic signature, since their velocity vectors should 
point almost directly away from the Galactic center (apart from small changes induced by 
the non-spherical component of the Galactic potential). 

Hypervelocity stars produced at the Galactic center should also be distinguishable from 
runaway stars produced by supernova explosions in close binary systems or by close encoun- 
ters in star clusters (Leonard 1993), which will have quite different velocity distributions. 

At a distance of lOkpc, a solar-type star traveling at lO^kms"-*^ has apparent magnitude 
V = 19.8 and proper motion 20masyr^^. Ground-based surveys can achieve proper-motion 
accuracies of ~ Imasyr^^ with multiple images over baselines of a few years (Alard 2003). 
Planned ground-based surveys such as Pan-STARRS and the Large Synoptic Survey Tele- 
scope, with limiting magnitude V ~ 24, could easily detect and measure the proper motion 
of solar-type hypervelocity stars over much of the Galaxy. The GAIA spacecraft will survey 
the whole sky to F = 20 with proper-motion accuracy of ~ O.lmasyr"^ or better at its 
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limiting magnitude. 

In Hills' mechanism, one component of the binary star is ejected, and the other is thrown 
onto a high-eccentricity orbit around the central BH (Gould & Quillen 2003). We expect that 
the number of young (10^ yr old) stars that are the remnants of tidally disrupted binaries 
should be ~ 100(77/0.1) at any time. 

Hypervelocity stars produced by Hills' mechanism should be approximately spherically 
distributed, while the spatial distribution of the hypervelocity stars due to interactions with 
the BBH should be flattened in the orbital plane of the BBH (see Zier & Biermann 2001). 

Support for this research was provided in part by NASA through grants from the Space 
Telescope Science Institute, which is operated by the Association of Universities for Research 
in Astronomy under NASA contract NAS5-26555. 
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